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Constraining the polarized gluon PDF in pp collisions at RHIC 

F. Ellinghaus, for the PHENIX and STAR Collaborations 
University of Colorado, Boulder, CO 80309, USA 

The main focus of the physics program at PHENIX and STAR that makes use of RHIC's polarized proton beams 
is to figure out how and if at all the gluons inside protons are polarized, or to put it another way, do the spin 1 
gluons prefer to have their spins aligned or anti-aligned with the spin of the proton, or do they just not care? This 
question is an important part of the more general question of how the constituents of protons, gluons and quarks, 
conspire to make up the overall spin 1/2 of the proton. Measurements of, e.g, jet and hadron, production cross-section 
differences between the two cases where the two polarized protons colliding have their spins aligned and anti-aligned 
are sensitive to the gluon polarization, which is encoded in the spin-dependent parton distribution function (PDF) 
for gluons, Ag{x). 

In the most naive model, the so-called constituent quark model, in which the protons consist of three static quarks, 
two of the three spin 1/2 quarks will align their spins with the spin of the proton, the third one is anti-aligned, 
and voila, the three quarks account for the spin 1/2 of the proton. However, after pioneering experiments at SLAC 
starting in the mid-1970's it was the European Muon Collaboration (EMC) at CERN in the late 1980's that found 
the contribution of the quark spins to the spin of the nucleon to be small and actually consistent with zero within 
the uncertainties [ij. The most recent results by the HERMES Collaboration at DESY [i] and the COMPASS 
Collaboration at CERN [sl show that the quarks contribute about 30% to the spin of the nucleon, leading to our 
current understanding that a substantial fraction of the proton's spin has to arise from the quark orbital angular 
momentum and/or the gluons. 

The basic principle of these measurements is to use deep inelastic scattering (DIS) of longitudinally polarized 
leptons off longitudinally polarized nucleons, which are usually embedded in a gaseous or solid state target. The 
polarized lepton will interact with a quark inside the polarized nucleon via the exchange of a spin 1 virtual photon. 
This photon can only be absorbed by the quark when the spins of the quark and the photon are anti-aligned due 
to angular momentum conservation. The cross section for this process can be measured for the two cases where the 
spin of the leptons and the nucleons is aligned and anti-aligned, and depending on which cross section turns out to 
be bigger, one can learn if the quark spins tend to be aligned or anti-aligned with the spin of the nucleon. This 
cross section difference is usually measured as a function of Q^, the negative four- momentum squared of the virtual 
photon, and of x, the fraction of the momentum of the (fast) proton carried by the struck parton, and it is given in 
terms of the so-called polarized structure function gi(x, Q^). 

The possible contribution of the gluon spin AG{Q^) = dxAg{x, Q^) to the nucleon spin has first been studied 
in DIS. The problem with studying gluons in DIS is that the electromagnetic probe, the photon, cannot directly 
couple to the gluon since gluons carry no electric charge. An indirect way of accessing the polarized gluon PDF 
Ag(x) is via next-to-leading order (NLO) fits to the above mentioned structure function gi{x, Q^). The information 
on the polarized gluon PDF is contained in the dependence of gi on at a given x, which is nearly constant. To be 
precise, it is the difference from the constant behavior, the so-called scaling violation, which contains the information. 
However, since all measurements of gi have been done at fixed-target experiments the range in covered is very 
limited, leading to rather large uncertainties for Ag(x). A more direct way to measure Ag{x) at the DIS experiments 
is via a NLO process called photon-gluon fusion (PGF), in which the virtual photon and the gluon interact via the 
creation of a quark-antiquark pair. However, the available results so far suffer from the fact that the energy scale at 
the fixed-target experiments is rather low. 

Thus, the motivation to study the gluon in pp collisions at the STAR and PHENIX experiments at RHIC (Rel- 
ativistic Heavy Ion Collider) at BNL is quite natural. The involved energies are high and the gluon participates in 
leading order since the basic processes of interest are the scattering of quarks off gluons and of gluons off gluons. 
Similar to the case of gi the cross section difference is measured between the two cases where the two polarized 



34*'' International Conference on High Energy Physics, Philadelphia, 2008 



0.1 
0.08 
0.06 
0.04 
0.02 


-0.02 



4 2006 STAR Preliminary 

— GRSV-std 

— GRSV Ag=g 
- . GRSV Ag=0 

— GRSV Ag=-g 
■-■ GS-C 




30 35 

(GeV/c) 



=„-j0.05 
< 














: PHENIX Preli 


minary 








0.04 
















: • Runs (9.4% scaling error) 








0.03 


■ Run6 (40% scaling error) 










0.02 










GRSV- 


std 


0.01 


L 1 






GRSV 


AG = 




-0 










GS-C 




-0.01 


Scaling errors not 


Included 








-0.02 


Online Polarization values 
used in Run6 








-0.03. 


-< , , , 1 < < , , 1 , < < < 1 , 


,1 


, , 1 


1 1 1 1 






2 3 4 


5 


6 


7 8 3 
(GeV/c) 



Figure 1: Double helicity asymmetry for mid-rapidity inclusive jet (left panel) production from the 2006 (preliminary, see 
Ref. for final results from 2005 data) data and tt" (right panel) production from the 2006 (preliminary") and 2005 data ^ 
at y/s = 200 GeV as a function of pr. The results for the asymmetries are compared to pQCD calculations \10^ based on 
different sets of polarized PDFs [ill. Il2|. See text for details. 



protons colliding have their spins aligned and anti-aligned. To be precise, the difference divided by the sum, the 
so-called double helicity asymmetry 
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is measured, where the cross section ct++ (ct^^ ) describes the reaction where both protons have the same (opposite) 
helicity. The spin-dependent term is given on the rhs of Eqn. [l] where A/a, A/{, represent the spin-dependent 
PDFs for quarks (u,d,s) and gluons, which, from the NLO fits to gi as mentioned above, are well known for quarks 
but not well known for gluons. The fragmentation functions represent the probability for a certain parton 
c to fragment into a certain hadron h, and thus they are not needed in the case of jet production. The spin- 
dependent hard scattering cross sections Act have to be calculated using perturbative QCD. Thus, the applicability 
of perturbative QCD in the kinematic regime of the A^l measurements is a crucial prerequisite. This has been 
checked by comparing the results from NLO pQCD calculations using the well known unpolarized PDFs for quarks 
and gluons with measured unpolarized cross sections for, e.g., the inclusive production of jets or neutral pions [5|. 
Good agreement has been found giving confidence in proceeding with measurements of the double helicity asymmetry 
(Eqn. [T]), which experimentally translates into 
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where A^++ [N^ ) and £++ [L^ ) are the experimental yield N and luminosity L for the case where the beams have 

the same (opposite) helicity. The relative luminosity R is just the ratio of the two luminosities. The polarizations 
of the two colliding beams at RHIC are denoted by Pb and Py- The degree of polarization is determined from the 
combined information of a pC polarimeter using an unpolarized ultra-thin carbon ribbon target, and from pp 
scattering, using a polarized atomic hydrogen gas-jet target The average polarization value for the data from 
2005 (2006) is 49% (57%). 



^Final results are now available in Ref. [^. 
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Figure 2: Confidence level of the STAR jet data (left panel) and profile for the PHENIX tv'^ data (right panel) resulting 
from a comparison to four GRS V scenarios as well as for several other values of AG at the input scale . 



The double helicity asymmetry as a function of the transverse momentum px for jet production at STAR and tt*^ 
production at PHENIX are shown in Fig. [TJ The results are compared to NLO pQCD calculations [10| , using two 
sets of polarized PDFs. In addition to the GS-C PDFs [lH, the GRSV PDFs [12] (GRSV-std) have been modified by 
assuming the polarized gluon PDF to be zero (AG — 0), equal (AG = G), or opposite (AG = — G) to the unpolarized 
gluon PDF at the input scale. It is apparent from Fig.[l]that the maximum and minimum GRSV scenarios, which are 
not even shown anymore on the PHENIX plot, are ruled out by both experiments. While not excluded, GRSV-std 
is disfavored by the data as can be seen in Fig. [21 where for the STAR jet data (PHENIX 7r° data) the confidence 
level (x^ profile) is shown for the four GRSV scenarios as well as for several other values of AG at the input scale. 
Note that no uncertainties have been calculated for the GRSV set of PDFs, and that they in general are rather large 
for the poorly constrained Ag{x) from the DIS data as discussed above (see, e.g., Ref. [13[), hence the RHIC results 
remain consistent with the results from DIS. 

Based on Fig. [2] a small or maybe negative contribution of the gluon spin AG to the nucleon spin appears to be 
most likely; however, this is a model dependent result. The present inclusive measurements cannot determine x but 
are sensitive to the gluon contribution integrated over the range 0.05 < a; < 0.2. Thus, in case the gluon contribution 
is large outside the presently accessible a;~range it could still account for all the "missing" spin, even though the 
contribution in the limited accessible x range is rather small. An example is the predicted asymmetry based on the 
GS-C PDF also shown in both panels of Fig. [1] It is consistent with the present data, but the GS-C PDF gives a 
large contribution of gluons to the nucleon spin due to the fact that it is large and positive at small x and becomes 
negative above x « 0.1. In other words, the gluon contribution is small in the measured region since it has a node, 
but the positive contribution at smaller x could account for all the spin due to the abundance of gluons at small x. 

These features can also be seen in Fig. |31 which shows a recent NLO pQCD fit (DSSV [l4|), including for the 
first time the pp data in addition to DIS and semi-inclusive DIS in the extraction of the polarized PDFs, a so-called 
global fit. In the x range accessible at RHIC the DSSV fit yields smaller values for xAg{x) than the earlier GRSV 



or DNS 15| fits with much reduced uncertainties. However, at smaller values of x the uncertainties are large and 
thus allow for a large xAg{x). A lower (higher) range in x can be probed runn ing RHIC at higher (lower) energies. 
Results on tt" production at PHENIX at ^/s ~ 62.4 GeV are already available [16|, while running at \/s = 500 GeV 
is foreseen to start in 2009. In addition, measurements at forward rapidity will allow access to a lower range in x. 
Also, with more data available in the future it is possible to advance from inclusive measurements, in which only one 
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Figure 3: Polarized gluon PDF from a fit [f^l that included for the first time the pp data (in addition to DIS and semi-inclusive 
DIS), compared to the earlier fits GRSV and DNS [li]. Alternative DSSV fits corresponding to Ax^ = 1 (Ax^/x^ = 2%) are 
given by the green (yellow) bands. Figure taken from Ref. [f^]. 



particle in the final state is detected, to measurements in which two or more leading particles or jets are detected 
in the final state, allowing for a much better handle on the ^-dependence. Since in DIS the momentum fraction 
X can be measured directly, a precise measurement of the functional from of the polarized gluon PDF At; (a;) will 
be possible at a proposed polarized electron-proton collider called EIC [l7|, which would be able to overcome the 
difficulties arising from the low energy scale and limited range at fixed-target experiments. 
This work is supported in part by the US Department of Energy. 
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